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2 I We have presented a complete description of classical dynamics generated by the 

O ' Hamiltonian of quadrupole nuclear oscillations and identified those peculiarities of quan- 

^ ■ tum dynamics that can be interpreted as quantum manifestations of classical stochas- 

ticity. Particular attention has been given to investigation of classical dynamics in the 
O I potential energy surface with a few local minima. A new technique is suggested for 

^ I determination of the critical energy of the transition to chaos. It is simpler than criteria 

of transition to chaos connected with one or another version of overlap resonances cri- 



terion. We have numerically demonstrated that for potential with a localized unstable 
region motion becomes regular at the high energy again, i.e. transition regularity-chaos- 
', regularity {R — C — R) takes place for these potentials. The variations of statistical 

properties of energy spectrum in the process of R — C — R transition have been studied 
in detail. We proved that the type of the classical motion is correlated with the structure 
of the eigenfunctions of highly excited states in the R — C — R transition. Shell structure 
destruction induced by the increase of nonintegrable perturbation was analyzed. 
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1 Introduction 

1.1 Formulation of the problem 

Quite often in the different divisions of physics statistical properties are introduced by 
means of postulates or hypotheses that are self-justified within some limiting case. Such 
an approach suggests an existence of a priori mechanism of randomness the nature of 
which lies usually outside the theory under discussion. For systems with a small number 
of degrees of freedom the results obtained in the thermodynamic limit look at least 
questionable. At the moment, however, one can consider as rigorously established fact 
an existence of such dynamical systems with a small number of degrees of freedom 
for which under certain conditions classical motion is distinguishable from random no 
matter what the definition of this notion we are using [|l], ^ ^ H, ^, |[. Typical features 
of these systems are nonlinearity and absence of both external source of randomness 
and dissipation the mechanism of which is equivalent to the existence of random forces 
on a molecular level. Thus, using such synonjnns for the term "random" as "chaotic", 
"stochastic", "irregular" one can state that there exist such systems with a finite number 
of degrees of freedom, for which these notions express adequately internal fundamental 
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properties that comprise an important and interesting subject for investigation. For the 
last 30-40 years has been taken place a rather difficult recognition of the fact that chaotic 
motion is as usual in systems with more than one degree of freedom as regular motion. 
Instances of a chaotic motion have been found practically in every field of physics, and 
their number permanently continues to grow [|l], |], ||, H, H, If. 

Recent progress in understanding of the chaotic aspects of nonlinear dynamical sys- 
tems is causing a rebirth of interest in the nuclear dynamics. Indeed, because of the 
richness of experimental data and sufficiently precision of the theory, the nuclear dy- 
namics provides a very useful realistic model for studying classical dynamical chaos and 
quantum manifestations of the classical stochasticity. 

Conception of chaos has been introduced in the nuclear theory within the last fifteen 
years 0, ||, |, |TU|, 0, 0, |T|, ^ |g, [T^, 0, 0]. This conception brought birth to 
the new notion in the nuclear structure |0, |ll|, nuclear reactions p|, [l^, [20[] , 



could resolve a sequence of the very old contradictions in the nuclear theory [|r^, A 
radically new universal approach to the problem of statistical properties of the energy 
spectra was developed on the basis of the general nonlinear theory of dynamical systems 
||23|| . Considerable advances have been made in the area of concrete nuclear effects: 
single particle in a deformed potential ||ll|], nuclear fission 0, Ericson's fluctuations ||24|| , 



dynamics of the 3a - linear chain , transition order-chaos-order in the roto- vibrational 
nuclear model ||53|, Q etc. Finally, straightforward observations of the chaotic regimes 
in the course of mathematical simulations of the heavy-ions reactions ^ are evidence in 
favour the general considerations. 

According to Baranger |TB[ there are two possible philosophies in nuclear physics. 

Philosophy I. Nuclei are complicated, and chaos comes out of this complication. We 
expect to find chaos almost everywhere in nuclear physics. The interesting information 
is contained in those few places in which chaos is absent. We must look for non-chaos. 

Philosophy II. Chaos is a property of simple systems; otherwise it's no fun at all. The 
interesting new information is to be found in those simple areas of nuclear physics which 
we used to think we could understand, but which turn out to be chaotic. We must look 
for chaos. 

The basis of present report is the philosophy of simple chaos - Philosophy II. Within 
the limits of this philosophy a general investigation of any nonlinear dynamical system 
involves the following steps. 

1. An investigation of classical phase space with the aim of detection of chaotic 
regimes; numerical investigation of the classical equations of motion . 

2. Analytical estimation of the critical energy for the onset of global stochasticity. 

3. Test for quantum manifestations of classical stochasticity (QMCS) in the energy 
spectra, eigenf unctions and wave packet dynamics. 
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4. A consideration of the interrelationship between stochastic dynamics and concrete 
physical effects. 

The basic subject matter of proposing report is to realize the outlined program, at 
least partially, for the large-amplitude quadrupole oscillations of nuclei (QON). We shall 
organize this paper as follows. 

In the first part (section 2.1-2.7) we present complete description of classical dynamics 
generated by the Hamiltonian of quadrupole oscillations. This description contains 

1. The scale analysis of the Hamiltonian (QON). 

2. The discussion of the topological peculiarities of deformation potential. 

3. The analytical estimations of critical energy of the transition to chaos. 

4. The construction of the approximate integrals of motion with the help of normal 
forms. 

5. The analysis of dynamics in the region of parameters space where the deformation 
potential has some local minima; the introduction of the concept of mixed state. 

6. The general conception of reconstruction of regular motion at high energies for the 
systems with localized region of instability; the transition regularity-chaos-regularity as 
a particular realization of this conception. 

7. Demonstration of the chaotic regimes in heavy-ions reactions. 

In the second part peculiarities of quantum dynamics of considered system, which 
can be interpret as the QMCS. For this purpose in this part of the report were discussed 

1. The Birkhoff's-Gustavson's quantum normal forms. 

2. The change of statistical properties of the energy spectra in R — C — R transition. 

3. The change of the structure of wave functions in R — C — R transition. 

4. The destruction and the reconstruction of the shell structure in R — C — R 
transition. 

5. The dynamics of wave packets. 

Thus, the distinctive feature of the suggested report is that all complex of questions 
connected with the investigation both classical stochasticity and QMCS is considered in 
the context of the unique dynamical system. 

As for circle of readers whom this report is counted on, then on the one hand we would 
like to pay attention of experts in all areas of chaos to perspective field of application 
of general theory of nonlinear dynamical systems, and on the other hand to persuade 
nuclear experts that ideology of simple chaos can be equally useful just as traditional 
statistical approaches. 

For convenience of the latter a brief summary of definitions and concepts used in 
description of stochastic motion in Hamiltonian systems will be given in section 1.2. 
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1.2 Classical stochasticity: some definitions 

The classical dynamics of Hamiltonian systems with N degrees of freedom is described 
by the canonical equations of motion 

rlT-f f)}-f 

P. = -^, g. = ^ {l^^^N) (1) 
oqi dpi 

where H {p, q) is the Hamiltonian of system. The function of dynamical variables F {p, q) 
is such that 



{F,H} = Q (2) 

where {., .} are Poisson's brackets and it is called the integral of motion. If there are 
independent integrals (1 ^ i ^ A^), such that {Fi, Fj} = 0, then the system (|l|) is 
integrable. 

For integrable system one can introduce such canonical variables action-angle (Jj, 6i) 
that Hamiltonian will be the function of only variables of action 

H{p,q)=H{I) (3) 

The region labeled by fixed / , to which the motion is confined, is therefore an N 
-dimensional torus in the 2N -dimensional phase space. Any phase point that lies on a 
given torus at any time remains on it for all future times, so the torus itself is invariant 
under the Hamiltonian flow, and is known as an invariant torus. The motion on an 
invariant torus is given by 



Oi (t) = u,t + Oi (0) (4) 
where angular frequency vector uoj can be written as 

In this case the motion is a quasiperiodic function of time with at least m independent 
frequencies. If the are not rationally related, that is, there are no integers s-i , such 
that 

^ SiU, = (s ^ 0) (6) 
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then the phase point passes arbitrary close to every point of the torus and it is not 
difficult to show that the time average of function is equal to the average of the function 
over the angle variables. Thus the motion on the torus is ergodic. If is satisfied for 
some nonzero integer vector s, than is a resonance, and N -tori are made up of tori of 
lower dimension, to which the motion is confined. Such -tori are degenerated. The 
particular case of integrable system is the system with separable variables. 

In the general case the system with two or more degrees of freedom is not integrable 
and can perform both quasiperiodic (regular) and stochastic motion. The distinctive 
feature of stochastic motion is instability, exhibiting in exponential divergence of close 
trajectories. If x(t) and x'(t) are two trajectories close at t = in phase space, then 
during stochastic motion 

A (t) = \x (t) -x'{t)\^A (0) e-^'W (7) 

at sufficiently small A (0). Quantity is called the maximal Liapunov exponent and is 
determined as 

ai{x)= lim filn^-^Vx = x(0)) (8) 



t^oo,A{0)^0 \ t A (0) 

The system with D -dimensional phase space is characterized by the set of A Lyapunov 
exponents ai , which are numbered in order of decreasing. The points belonging to 
trajectory possess equal values (Xj (x) . The motion is called stochastic if for trajectory 
CTj > and regular if cij = . 

Kolmogorov's entropy h is connected with Lyapunov exponents. For given stochastic 
trajectory 

h=J2^^ (9) 

At stochastic motion the correlation function 

B if, g, r) = {f{t + r)g (t)) - (/ (t)) {g {t)) (10) 

(corner brackets (...) denote the averaging on time, and f if) = f {p, q) and g{t) = g {p, q) 
tends to zero with the growth of t 

limS(/,(7;r) = (11) 

r— >oo 

This property is called the mixing. 
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Power spectrum S{uj) of the dynamic value / (t) is determined by the expression 



oo 



^(/;^) = ^ / B{fJ;T)e-'-^dT (12) 

— oo 

For the motion with mixing the power spectrum S{uj) is continuous, and for the regular 
motion it is discrete: 

S{f,uj) = Y,AkS{uJ-uJk) (13) 

k 

here form the countable sequence, and frequencies uj^ are the combinations of fre- 
quencies ujj of quasiperiodic motion. 

2 Regular and chaotic dynamics of 
the quadrupole oscillations 

2.1 Hamiltonian 



It can be shown ||26| that using only the transformation properties of the interaction, the 



deformational potential of surface quadrupole oscillations of nuclei takes on the form 

U (ao, as) = J2 Crnn {al + 2a^)™ (6a^ - ag)" (14) 



m,n 



where and 02 are internal coordinates of nuclear surface at quadrupole oscillations 
R {6, ^) =Ro{l + ao>^2,o {0, f) + ^2 [>^2,2 (^, ^) + 1^2,-2 (^, v)]} (15) 

Constants Cmn can be considered as phenomenological parameters, which within the 
limits of the particular models or approximations (for instance, ATDHF theory) can be 
directly related to the effective interaction of the nucleons in nucleus [^. Considering 



that at the construction of (n) only transformation properties of interaction are used, 
then this expression describes potential energy of quadrupole oscillations of a liquid 
drop of any nature, including only the specific character of a drop in the coefficients of 
expansion Cmn ■ 

Restricting ourselves with the members of the fourth degree in the deformation, and 
assuming the equality of mass parameters for two independent directions, we get the 
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following C^,y - symmetric Hamiltonian 

pI + pI 



H = ^-^!—^ + U{x,y;aAc), (16) 
2m 



d / o o\ ,/o 1q\ ^2 2\2 



where 

U {x, y- a, c) = ^ {x^ + y'^)+h \yx^y - ) + c (x^ + y^) 

a; = A/2a2, y = clq 

2 = Cio, & = 3Coi, c = C2o (17) 

The potential t/ (x, y] a, b, c) is the generalization of the well-known Henon - Heiles 
potential ||2^ with the one important difference that the motion in the potential ( |T7| ) is 



finite for all energies. This is particularly important for the quantum case, where the 
potential (0) ensures the existence of stationary states (instead of quasistationary ones 
as with the Henon- Heiles potential). 

Csi, - symmetry of potential surface becomes obvious in polar coordinates 

X = /3sin7, ?/ = /5cos7 (18) 

where /3 is the so-called parameter of deformation of axial symmetric nucleus, and 7 is 
the parameter of nonaxiality. 
In these coordinates 

U (/?, 7; a, b, c) = -a/?2 - -bf3^ cos 87 + c(3^ (19) 
2 3 

There are three possibilities of introduction of the typical unit of length /qj (^ = 1,2,3) 
for the Hamiltonian (|T6]) 

1) as the distance /qi, at which the contributions from harmonic and cubic terms 
become comparable 

2) as the distance /02 , at which the contributions from harmonic and biquadratic 
terms become comparable 

3) as the distance /qs , at which the contributions from cubic and biquadratic terms 
become comparable. 

Scaling of the principal physical values (scaling of the dynamic variables and energy) 
(x, y) = loi (x, y) , {p^,Py) = poi (Px, Py) ; E = eoiE (20) 
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for these three cases is determined by the following parameters 



^01 = 


a 
V 


P01 = \ 


/ a3 


^01 = 




lo2 = 


V 


-, P02 = 


/ a2 


^02 = 


c 


hs = 


h 

c 


P03 = 


/ 


^03 = 


C3 



(21) 



The reduced Hamiltonian for these three variants of scaling is the following 

HJp,,Py,x,y-W^-\ =^L±i + Udx,y,W) (22) 
\ ac J 2 

where 

Ui {x, y;W) = ^ {x^ + f) + (^^y - + ^ + ^')' 

f/2 (X, mW) = ]^ (X^ + f)+^ (^-^y - + (a;2 + ^2)2 

f/3 (x, W^) = ^ (x^ + f) + (^x^y - + (x^ + y2)2 (23) 

The transition i ^ k between the different variants of scaling is realized with the help 
of the substitution 

(x,y) {x,y) {px,Py) {Px,Pv) E ^ —E (24) 

hi POi ^Oi 

where the parameters of transformation are defined by the relations (pT]). 

At any variant of scaling the reduced Hamiltonian and corresponding equations of 
motion depend only on parameter W. This is the unique dimensionless parameter, which 
can be constructed from the dimensional values a, b, c. It means that for each trajectory 
of initial "physical" Hamiltonian (^) with energy E , corresponds the unique trajectory 
of one-parameter Hamiltonian (^) with energy E = E / e^i . While for each trajectory of 
the reduced Hamiltonian HiiW) with energy E corresponds the whole set of "physical" 
trajectories with energy E = SQiE, which are generated by Hamiltonians H{a, b, c) with 
parameters satisfying the conditions b"^ /ac = W. 
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Figure 1: The lines of the level of PES for different structurally stable regions. 
2.2 Geometry of potential energy surface 

Now, let's investigate the geometry of two-dimensional set of the potential functions 
U{x, y; W) . The set of solutions of the system of equations 

[/; = 0, Uy^ 0, det ^ = (25) 
where S is the matrix of stability 

(26) 

serves by separatrix in the space of the parameters and divides it into the regions, where 
the potential function is structurally stable. The number and the nature of critical points 
change at the transition through separatrix W — 16 and under the change of signW 
(i.e. under the change of sign a; c is always positive). The critical points of PES for 
each structurally stable regions are given in Table 2.2.1. 



Table 2.2.1. Number of critical points in different ranges. 





range 


critical points 


saddles 


minima 


maxima 


I 


< < 16 


1 





1 





II 


> 16 


7 


3 


4 





III 


w <0 


7 


3 


3 


1 



The corresponding lines of the level of the potential are represented in Fig.l. 



S 
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The coordinates of critical points are determined by 



xi = 0, yi = 




(27) 



X4-7 = ±V32/4-7, 2/4-7 = ( 1 ± - ^ 



where Oi = W,a2 = VW, 03 = 1. 

The region of the space of parameters < < 16 includes the potentials possessing 
the unique extremum - the minimum in the origin of coordinates, corresponding to the 
spherically symmetric equilibrium shape of nucleus. The region W > IQ includes the 
potential surfaces with four minima. The central minimum, with x = y = , corresponds 
to the spherically symmetric equilibrium shape of nucleus, but three peripheral minima 
correspond to deformed shape. At last, in the region IV < (a < 0) we run into the 
potentials, describing the nuclei, which are deformed in the ground state and even don't 
have the quasistable spherically symmetric excited state. This region can be divided 
into two subregions. In the space of parameters the boundary between these subregions 
W = —1/2 represents the geometric locus, where the both eigenvalues of matrix of 



stability turn into zero simultaneously. Seiwert,Ramayya and Maruhn |29| restored the 
parameters of the Hamiltonian of the quadrupole oscillations, including the sixth degree 
members in deformation for isotopes /^7-74,76,78,80 rjj^g experimental values of energy 
of the first 2+ - states for nuclei Tfr^"^'^^ indicate spherical shape of nucleus surface, while 
the probabilities of the electromagnetic transitions 2"*" —>■ O"*" and very low energies of the 
first rotational states (recalculating to one and the same number of nucleons, this energy 
for isotopes Kr'^'^ is essentially lower than the lowest from the known energy for nucleus 
p^24o _ ^2.8 kev ) imply the possibility of " superdeformation" ||3^. The nonlinear effects, 
which are connected with the geometry of PES, must be exhibited in the superdeformed 
nuclei even at the relatively low energies of excitation. The potential surfaces of Krypton 
isotopes, which are calculated in are presented in Fig.2. As it is seen, the inclusion 



into the consideration of the highest members of expansion in deformation leads to a 
considerable complication of the geometry of PES: for all considered Krypton isotopes 
we run into the PES of complicated topology with many local minima. It's clear, that 
in any degree in deformation of the PES C^u - symmetry is preserved. 
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Figure 2: The surfaces of potential energy of Krypton isotopes. 
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2.3 Critical energy of transition to chaos 



If we understand under the stochastization the process of appearance in the system of 
statistical properties in consequence of local instability, we obtain a tempting possibility 
of identifying the values of parameters, that lead to the local instability in the system, 
with the boundary of transition to chaos. Unfortunately, the criteria of the stochasticity 
of the similar type, i.e. based on the investigation of the local instability, posses the 
innate lack: the lost of stabihty of the regular motion does not obhgatory leads to chaos. 
Generally speaking, instead of this the transition to another more complicated type of 
regular motion is possible. Besides that, the statement, that the local instabilities define 
the global dynamics of the system is disputed. Separate details of the derivation of 
the concrete criteria of the stochasticity provoke several objections. In spite of these 
serious lacks the available experience allows us to state, that the criteria of the similar 
type in the aggregate with numerical experiment essentially facilitate the analysis of the 
many-dimensional nonlinear motion. 

A large variety of criteria of the stochasticity is based on the direct evaluation of 
the rate of divergence of the initially close trajectories {qi{t) , pi{t)} and {(hit) , P2{t)} . 
Linearized equation of motion for a divergence 



at) 



qi it) - q2 it) , V{t) 



Pi it) -nit) 



(28) 



assume the following form 




The stability of motion of the dynamic system, described by the Hamiltonian 




(31) 



is determined in N- dimensional case by 2N x 2A'" matrix 



r = 



o I 

Sit) 6 



(32) 
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where O and J is a zero and a unitA^ x matrices. If even one of the eigenvalues Aj of 
the matrix F is real, then the divergence of the trajectory increases exponentially and 
the motion is unstable. The imaginary eigenvalues correspond to the stable motion. The 
eigenvalues, and consequently the character of motion change with the time. 

The problem of the investigation of the stable motion can be essentially simplified 
pll , if we assume the possibility of a replacement time-dependent point of phase space 
'qiit) by time-independent coordinate q . It reduces the equations for variations (|29|) to 
a system of autonomous linear differential equations 

f=r/ ; f^=-S{l^U (33) 
An equation for Lyapunov exponents 1 which determine the character of motion 



det r - A/ 







(34) 



for the system with two degrees of freedom (the system of our interest) has the following 
solution 



Ai 



2,3,4 



± 



-(3 ± _ 4^ 



where 



/? = Sp'S = + Uyy, 

7 = det S 



TJ TJ _ r/2 



dx 



dy 



2 ^^xy 



dxdy 



(35) 



(36) 



Here we will assume that /? > 0. Then, providing that 7 > , Lyapunov exponents 
are purely imaginary and the motion is stable. With 7 < the pair of roots becomes 
real and it leads to exponential divergence of close trajectories, i.e. to the instability of 
motion. 



Now let's remind several known facts from the theory of surfaces |]32|. Gaussian 
curvature of a surface is equal to the ratio of the determinants of the second h and the 
first g quadratic forms 



K 



&11&22 - h 
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det 6 



9ii922 - 912 det g 
In particular, if the surface is given in the form of the graph z = U {x,y) , then 



(37) 



det 6 
det g 



TJ TJ _ r/2 

^xx^yy ^xy 



(38) 
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and therefore 

The Gaussian curvature K can be represented in the form of the product of the so-called 
principal curvatures di and d2 , which are the eigenvalues of the pair of the matrices h 
and g and are the solutions of the equation 

det {b -6g^ =0 (40) 

The sum of di + d2 is called an average curvature of the surface. We will enlarge on 
the geometric sense of the Gaussian curvature. Let's choose the orthonormalized frame 
{x,y,z) for given point XQ,yQ of surface, where z axis is a normal to the surface. Then 
locally the surface will be written as z = U{x,y) and f/^ = f/^ = in this point. 
Therefore 

S1S2 =K= UxxUyy - U% (41) 

We will consider three possible cases 

1) > 0, 5i > 0, (^2 > (minimum of the function U (x, y) with x = xo,y = yo , 

2) K > 0,61 < 0,62 < (maximum of the function U{x, y) with x = x^^y = yQ , 

3) K < 0, (5i < 0, (^2 > or vice versa (saddle point of the function U{x,y) with 
X = xo,y = yo . 

At i^' > the surface is locally situated on one side of tangential plane to the 
investigating point. At i^' < the surface necessarily crosses the tangential plane as close 
as possible to the point of tangency. If the Gaussian curvature is positive everywhere, 
then this surface is strictly convex. 

Recently, in different sections of physics a definite interest has grown to the surfaces 
which everywhere possess a negative curvature . Such surfaces in the neighbourhood 



of any point behave as in the neighbourhood of the hyperbolic singular point. Now let's 
return to the expression of Liapunov exponents in the case of two-dimensional potential 
surfaces. Comparing the (BB) and (|55|) expressions we notice that the g sign coincides 



with the sign of Gaussian curvature of the PES. This association suggests [3|, 35 1 the 
possibility of the existence of the following scenario of the transition from regular to 
chaotic motion, based on the investigation of Gaussian curvature of the PES. 

At low energies the motion near the minimum of the potential energy, where the 
curvature is obviously positive, is periodic or quasiperiodic in character and is separated 
from the instability region by the zero curvature line. As the energy grows, the "particle" 
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will stay for some time in the negative- curvature region of the PES where initially close 
trajectories exponentially diverge. At large time these results in the motion which imitate 
a random one and is usually called stochastic. According to this stochastization scenario, 
the critical energy of the transition to chaos, Ecr , coincides with the lowest energy on 
the zero curvature line 

Ecr = U^UK = 0) (42) 

In the subsequent discussion we will reference to this statement, as to the negative 
curvature criterion (NCC) . 

Now we shall concentrate our attention on one of the first (1959 year) widely used 
criterion of transition to chaos, the so-called overlap resonance's criterion (ORG) |^ (the 



criterion of Chirikov). The essence of this criterion is easier to explain by the example 
of one-dimensional Hamiltonian system, which is subjected to monochromatic periodic 
perturbation. This one is the simplest Hamiltonian system which assumes the chaotic 
behaviour 

H = Ho (p, x) + Fx cos nt (43) 

For unperturbed system we can always introduce the variables action-angle (/, 6) in 
which 

oo 

H = Ho{I)+ Xk{I)cos{ke -Qt) (44) 

fc=— oo 

where 

^k{I) = ^ [ dee'^'x{I,9) (45) 
ivr J 



In new variables the scenario of stochasticity, on which the overlap resonance's cri- 
terion (ORG) is based on, is the following. An external field which is periodic in time 
induces a dense set of resonances in the phase space of a nonlinear conservative Hamil- 
tonian system. The positions of these resonances, Ik , are determined by the resonance 
condition between the eigenfrequency uju (I) = and frequency of the external pertur- 
bation, Q . For very weak external fields the principle resonance zones remain isolated. 
As the amplitude F of external field is raised, the widths Wk of the resonance zones 
increase 



Fxk \ 



u;'{I) 



(46) 

i=h 
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and at F > Fcr resonances overlap. When this overlap occurs, i.e. under the condition 

^{Wk + Wk+i) = \Ik-Ik+i\ (47) 

it is said that there is transition to a global stochastic behaviour in the corresponding 
region of the phase space. In other words the ORG postulates, that the last invariant 
KAM surface, which separates the neighbouring resonances, is destroyed in the moment 
of contact of unperturbed separatrices of these resonances. In other words , the averaged 
motion of the system in the neighbourhood of the nonlinear isolated resonance on the 
plane of the variables action-angle is similar to the particle behavior in the potential well. 
Several isolated resonances correspond to several isolated potential wells. The overlap 
of the resonances means, that there is such an approach of the potential wells, wherein 
the random walk of a particle between these wells is possible. 

The outlined scenario can easily be "corrected" for the description of the transition 
to chaos in the conservative system with several degrees of freedom. The condition of 
the resonance between the eigenfrequency and the frequency of external field must be 
replaced by the condition of the resonance between the frequencies, which correspond to 
different degrees of freedom 

= (48) 

The role of the amplitude of the external field in this case plays the intensity of the 
interaction between different degrees of freedom , i.e. the measure of nonlinearity of the 
original Hamiltonian. This parameter is usually the energy of system. 

This method must be slightly modified for the systems with the unique resonance. In 



this case the origin of the large-scale stochasticity is connected [0 with the destruction 
of the stochastic layer near the separatrix of this unique resonance. The essence of the 
modification consists in the approximate reduction of the original Hamiltonian in the 
neighbourhood of resonance to the Hamiltonian of nonlinear pendulum, which interacts 
with periodic perturbation 

H {v, X, r) = -f ^ — M cos x — P cos k{x — t) (49) 



The width w of the stochastic layer of the resonance is equal to[pS 



w^P^ ^ Mp^^+^ (50) 
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where 



If p/M has the power , then 



where 



\ = 2k + l-s. 



At 



Pi 



A 



1-(1 + A) 



-12 



(51) 

(52) 
(53) 

(54) 



the function w{r) has a point of inflection. The fast growth of w allows us to determine 
the thresholds of the destruction of the stochastic layer as a value 

n 2 



ps = A" 



;i + A) 



1 



■/2 



(55) 



This value p is such that the tangent to the point of the function w crosses the axis p 

All three described criteria of stochasticity will be used below to determine the critical 
parameters of the transition to chaos. 



2.4 Numerical results versus analytical estimations 

Now let us turn to the analysis of the solutions of the equations of motion, which are 
generated by the Hamiltonian (|16]). As mentioned above in the section 2.2, the geometry 
of the PES for the regions 
10<W <16 

II W> 16 

III W <0{a< 0) 

is essentially different. Undoubtedly the specific character of the PES must be mani- 
fested in behaviour of the solutions of the equations of motion. Therefore we shall analyze 
each of the mentioned regions separately. Notice, that parameters of the Hamiltonian 
a,b,c were estimated in different phenomenological models [p9|. They change in such 

wide limits that for real nuclei parameter W = ^ /^^^ can belong to each of the regions 
I,II,III. 
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2.4.1 Region < < 16 : potentials with unique extremum 

It is the simplest region for the analysis in which good agreement between the boundary 
of the transition to chaos, observed in numerical experiments and analytical evaluations, 
can be obtained with the help of the simplest of the criterion of the transition to chaos - 
the negative curvature criterion. Once more we remind, that according to this criterion 
the mechanism of the generation of the local instability consists in hit of the particle in 
the region of negative Gaussian curvature of the PES and critical energy of the transition 
to chaos coincides with the minimum value of the energy on the zero Gaussian curvature 
line. The equation of the last one in re-scale coordinates (further, we shall use the third 
variant of scaling (P3|)) is 



1 W -4 



{x^ + f) + (x' + fy- (^x^y - 



and potential energy on the zero curvature line is 



S']=0 (56) 



= Ssff^ - ^ + f)+2{S:^ + sr (57) 

As it is seen, the zero curvature line conserves the symmetry of the PES (see Fig. 3.). 
Therefore, the minimum of the energy on the zero curvature line must lie either on 
the straight line x = or on the straight lines, obtained from it with the help of the 
transformations of the symmetry of the discrete group Csi^. 

At < < 4 Gaussian curvature of the PES is positive everywhere. At 4 < ly < 12 
the subregion of the negative curvature is localized on the straight line (in the plane) 
X = in the interval 



<y<-- 1- Jl-— (5J 





and in the subregion 12 < ly < 16 the negative curvature appears at ?/ > in the 
interval 

1 / I 1 O \ 1 / I 1 o \ 

(59) 

In Fig. 4. the profiles f/ (x = 0, y) of the PES are shown for the three considered 
subregions, and the intervals of the negative curvature are shaded. 




19 




Figure 3: Isolines of the PES for = 13 . Zero curvature line K (x, = (dashed 
hue) is shown on the background of the level hues, the range K (x, y) < is shaded. 

Thus, according to the considered scenario of the stochastization in the subregion 



< < 4 the motion must remain regular for all energies. According to the (P2| ) in 
the neighbourhood of the energy 



Ecr = f/min (K = 0) = U (x = 0, , T/erl = U " " 777 (^O) 




the transition to the global stochasticity should be expected in the subregion 4 < W < 
12. Finally, in the sub region 12 < < 16 this transition must be observed in the 
neighbourhood of the energy 



Ecr = f/min {K = 0) = U (x = 0,ycr2),ycr2 = ^2 ~ \l ~ W ) ^^^^ 

In this case we have used that for all 12 < < 16 f/(a; = 0, ?/cri) > U{x = 0, ycr2) ■ 

These predictions must be compared with the numerical solutions of the equations 
of motion generated by the Hamiltonian (|2^). 

Perhaps, the simplest numerical method of the detection of stochasticity is the anal- 
ysis of the Poincare surfaces of section. As it is well known [40] the Hamiltonian system 
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Figure 4: The profiles U {x = 0, y) of tlie PES for tlie sub regions 0<iy<4, A <W < 
12 and 12 < 1^ < 16 . The regions of the negative curvature are shaded. In the top 
part of the Figure - the Poincare surfaces of section for the energy values indicated by 
arrows, are represented. 

with N degrees of freedom is integrable only in the case of the existence of N indepen- 
dent univalent integrals of motion. If the number of the integrals of motion is less than 
the number of degrees of freedom, then the dynamic chaos is possible in the system. In 
the case of the existence of one or more additional (except energy) integrals of motion, 
the points of intersection of the phase trajectory with the arbitrary chosen plane (it is 
just Poincare surface of section) lie on the variety, the dimension of which is smaller than 
2N — 2 . Otherwise they will fill all the 2N — 2 dimensional isocncrgctic volume (all 
2 A'" — 2 isoenergetic surface). Therefore the analysis of the Poincare surfaces of section 
allows us to establish the fact of the existence of the additional integrals of motion and, 
consequently, to elucidate which type of motion is realized in the system with the given 
initial condition. The analysis of Poincare surfaces of section is especially effective for 
the system with two degrees of freedom, the phase space of which is four-dimensional. In 
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view of the conservation of energy, the trajectory of the particle hes on three-dimensional 
surface H {px,Py, x, y) = const . Excepted one of the variables, for example , we shall 
consider the points of intersection of phase trajectory with the plane x = . In common 
case they will be chaotic distributed on some part of the plane (p^, y) restricted by separa- 
trix. In the case of existence of the additional integral of motion / {px,Py,x,y) = const 
the totality of the consecutive intersections with the chosen plane lies on some curve 
f{Py,y) = const . By contrast, chaotic trajectories are identified by the fact that they 
show no orderly pattern on the Poincare surface of section. As we move from a regular 
regime into the chaotic one, disorderly trajectories seem to appear first near separatrices 
between different types of regular orbits. They become visible only in small regions, but 
as we move further into the chaotic regime, they cover most or almost all of the surface 
of section until eventually no regular trajectories are visible. 

The analysis of Poincare surfaces of section, obtained by the numerical integration 
of Hamiltonian equations of motion for the PES with unique extremum {0 < W < 16) , 
leads to the following results: 

1) At low energies {E « Ecr) in the neighbourhood of minimum for all considered 
values W , the motion remains regular. The regularity of the motion with low energies 
is a straight consequence of the KAM theorem which states, that the majority 
of the regular trajectories of the unperturbed systems remain regular under sufficiently 
small perturbation. It is clear that this is accorded with the positiveness of Gaussian 
curvature in the neighbourhood of any minimum. 

2) In the interval < W < A for all energies the motion remains regular (see Fig. 
4a.). This can be explained by the positiveness of the Gaussian curvature of the PES in 
this region of the parameter W . 

3) In the interval 4 < < 16 as the energy increases, the gradual transition from 
the regular (quasi periodic) motion to chaotic one (see Fig. 4b.) is observed. Moreover 
in the subregion 4 < W < 12 the critical energy, observing in the Poincare surface of 
section, is close to (pD|); as in the subregion 12 < < 16 the transition to the global 
stochasticity is observed in the neighbourhood of the energy, calculated according (|6l|) . 
This effect at 12 < 14^ < 16 is connected with the beginning of the region of negative 
curvature at y > located below in the energy (see Fig. 4c.). 

It is necessary to make one important remark. Analysis of the Poincare surfaces of 
section allows to introduce the critical energy of the transition to chaos, having deter- 
mined it as the energy in which the part of phase space with chaotic motion exceeds 
certain arbitrary chosen value. Similar indetermination is connected with the absence 
of the sharp transition to chaos for any critical value of the perturbation to which inte- 
grated system is undertaken. Therefore a certain caution is required in comparing the 
"approximate" critical energy, obtained by any variant of numerical simulation, with 
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the "exact" value obtained with the help of analytical estimations, i.e. on the base of 
different criteria of stochasticity. 

Based on this remark we can say, that in the case of one-well potentials, the negative 
curvature criterion allows to make reliable predictions about the possibility of the exis- 
tence of chaotic regimes in the considered region of the parameters; and also to evaluate 
the region of energies at which the transition regularity-chaos is performed. 



2.4.2 Region W > IQ : the potentials with a few local minima 

Now we are proceeding to the analysis of numerical solutions of the equations of motion in 
region W > 16 . The geometry of the PES, which is more complicated in comparison with 
the potentials which have the unique extremum (see Table 2.2.1), assumes the existence 
of several energies of the transition to chaos even for the fixed set of parameters of the 



potential. It means, that for such potentials the so-called mixed states |^ must be 
observed: at one and the same energy in different minima the various dynamic regimes 
are realized. The Poincare surfaces of section at different energies for Hamiltonian (p^) 
with W = 18 are presented in Fig. 5. This value W provides equality of depths for the 
central and peripheral minima. Taking into account C^i, symmetry of the PES, only one 
peripheral minimum is presented in this Figure. The motion represented in Fig. 5. a has 
clearly defined quasiperiodic character both for the central (left minimum) and for the 
peripheral (right) minima. Special attention must be given to the distinction in structure 
of the Poincare surfaces of section for different minima: the complicated structure with 
several fixed points in left minimum and simple structure with the unique fixed elliptic 
point in right one. The gradual transition to chaos is observed with the increasing of 
energy, however the change of the character of motion of the trajectories, localized in 
certain minimum, is essentially different. Whereas there is the gradual transition to chaos 
for the left well even for the energy approximately equal to one-half of the saddle energy 
(Fig. 5.b), and for the energy equal to the saddle one (Fig. 5c), practically all initial 
conditions correspond to the chaotic trajectories, the motion remains quasi periodic in 
the (second) right minimum at the same energies. In this minimum the transition to 
global stochasticity takes place only at the neighbourhood of the saddle energy. In the 
right well the significant part of phase space corresponding to quasi periodic motion 
conserves even at the energy essentially exceeding the saddle one (Fig.Sd). Figure 6. 
shows a comparison of the critical energies of the transition to chaos obtained according 
to the negative curvature criterion and with the help of the analysis of the Poincare 
surfaces of section. On the base of this comparison we can do the following conclusions. 



1) The mixed state is observed for considered potential for region of energies l/2Es < 
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Figure 5: The Poincare surfaces of section for the motion in the potential with W — 18 
at different energies: at the energy less then critical, E <C Ecr , (b) at the energy above 
critical, E > Ecr , (c) at the energy equal to saddle, E = Es .and (d) at the energy 
above the saddle, E > Es ■■ The profile U {x = 0,y) of potential is presented below. 
The range of < is shaded. 
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W=12 W=16 
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Figure 6: The critical energy of the transition to chaos for region: the negative curvature 
criterion (rigid hne) and with the help of the analysis of the Poincare section. 



E < Es{Es— saddle energy). 

2) The critical energy, defined according to the negative curvature criterion for the 
left well of potential, is in a good agreement with obtained one by means of numerical 
integration of the equations of motion and contradicts to the situation for the right well, 
where the numerical simulation detects large-scale chaos only at attainment of the saddle 
energy. 

The mixed state, which is shown above for the potential of quadrupole oscillations, is 
the representative state for the wide class of two-dimensional potentials with a few local 
minima. We shall show this by the example of the polynomial potentials of the degree not 
higher than six, which are symmetric relatively to the plane x = . But even with such 
restrictions, the possible set of the potential forms, depending (generally speaking) on 
12 parameters, is too big. Remaining generality, we shall use the methods of the theory 
of catastrophe |H| in order to reduce amount of calculations. According to the last one, 
a rather wide class of polynomial potentials with the several local minima is covered by 
the germs of the lowest umbilical catastrophes, of type D^,D^,Di , subjected to the 
definite perturbations. Let us notice, that the potential of type Henon-Heiles coincides 



with the elliptic ombilic |Q with precision to linear terms of perturbation. 



The mixed state is observed for all considered potentials of umbilical catastrophes in 



25 



the interval of energies E^^ < E < Es (here Ecr is the critical energy of the transition 
to chaos determined by negative curvature criterion). The transition to chaos (contrary 
to the negative curvature criterion) is observed only at the reach of the saddle energy 
for the minima possessing unique fixed elliptic point in the Poincare surface of section, 
as in the case of the potential of quadrupole oscillations. This contradiction makes one 
to refer to the criteria of stochasticity described in the section 2.3., and based on the 
theory of non-linear resonance. 

Let us consider, for example, the Hamiltonian with the PES which represents the 
germ of the catastrophe -D5 with quadratic perturbation 



h = 1{pI + pI) + 



4 I 2 I 2 2 
-y +x y + ax -y 



(62) 



The geometry of the two-dimensional one-parameter potential U{x, y; a) at a > \/2 is 
determined by five critical points : by two minima of equal depth (which thereafter will 
be named left and right wells) and by three saddles. The energy of the saddle, situated 
at the origin of the coordinates and separating the wells, does not depend on a and 
is equal to zero. Therefore, the classical motion will be localized at the separate well 
at the negative energies. The energies in all the saddles coincide with each other at 
a = 2 (this case will be investigated in detail). Now let us estimate the critical energy 
of the transition to chaos by the negative curvature criterion. In this case the problem 
is reduced to the search of the conditional extremum - minimum of the potential energy 
on the zero curvature line. The last one is described by equation 

(y + a) (3y2 - 2) - 2x2 = (63) 

At a = 2 we come to the value of the critical energy, which is equal for the both wells 

E,r = f/min {K = 0) = 

The results of numerical integration of the equation of motion generated by the 



Hamiltonian ( p^ ) and presented in Fig. 7 qualitatively coincide with the results obtained 
for the Hamiltonian of quadrupole oscillations (Fig. 5). Well manifested mixed state is 
observed at the energies E ^ ^2^8 iii Poincare surfaces of section. 

For the explanation of this phenomenon we shall consider the dynamics in the vari- 
ables angle-action. The Hamiltonian (^) in the system of coordinates with the origin 
in the left (top sign) and in the right (bottom sign) has the following form 

H=^{x' + oofx^) + + ^W) + x^y T V2l/' + \y' (64) 
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Figure 7: The results of numerical integration of the equation of motion generated by the 
Hamiltonian with PES which represents the germ of the catastrophe with quadratic 
perturbation (|5^. In the top part - the profile U {x = 0,y) of potential. Below - the 
Poincare surfaces: (a) at the energy less then critical, E <^ Ecr , (b) at the energy 
sligthly above critical, E > Ecr and (c) at the energy equal to saddle, E = Es ■ 



where uji = [2 (a =f ^2)] ,uj2 = 2. 

Now we perform the canonical transformation to the variables action-angle of the 
oscillator part of the Hamiltonian 

^= {^) ^^cos(/.i, y= ^^cosy.2, ^g^^ 

X = {2IiuJi) ''2cos(y9i, y = {2I2UJ2) ^'^cos(p2- 



In these variables the Hamiltonian (^4]) has the following form 

H (h, I2, ^1, '^2) = Ho (Ii, I2) + fm^mAh, h) cos (niiLpi + m2ip2) 

mim2Gy (66) 

y: [0,1], [0,2], [0,3], [0,4], [2,1], [2,-1] 

where 

3 P 

Ho (/i, h) = huJi + I2UJ2 + ^4 - 1 

O UJn 
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1/, 



8 2uJi \uj2 J 

The term with indexes r = (ri,r2) is called the resonance term for the given value of 
energy E , if there are variables action (/{■, IJ^) such that E = Hq (/[, I2) '■ 

nuj,{IlJl)+r2UJ2{IlJl) = Q (68) 

where 

0Oi = 9H^/QI^,i = 1,2. (69) 

It is possible to avoid problems of "small dominators" to perform canonical transfor- 
mation to the new action-angle variables which eliminate the dependence of the angle 
in the lowest order in small parameter identified with the energy at the energies corre- 
sponding to the finite motion (— 1 < < 0) , if we are situated sufficiently far from the 
resonance, i.e. for all mi, 1712, 

\miLJi + m2LJ2\ ^ /mima (70) 

The result of this procedure is the overdetermination of an integrable part Hq{Ii, I2) 
of the initial Hamiltonian and the extension of the set J members depending on angles. 
In this situation it is possible to meet one of the following three cases: 

1) the resonance terms before, absent in the region of energies of our interest; 

2) the unique resonance term appears; 

3) several resonance terms appear. 

In the first case we perform a new canonical transformation and continue the pro- 
cedure until we meet the situation 1) or 2). In the second case the critical energy, at 
which the transition to large-scale stochasticity takes place, can be determined by the 
destruction stochastic layer method [^], while the overlap resonances criterion can be 
used for the determination of the critical energy in the third case. 
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Small detuning near the resonance 



\miUJi + m2UJ2\ ^ /, 



r?iim2 



(71) 



can be compensated by terms of the highest order obtained with the help of repeated 
canonical transformation of non-resonance terms. In the case when a = 2 for the reso- 
nance (2,-1) the condition (^) of small detuning is performed and procedure, mentioned 
above, leads to 



+ 



3^2 



4uji + 5 



3,2 

8 ^ 32cj2(cji + l) 



UJl {uJi + 1) 



(72) 



In the left well among terms depending on angles we leave only two terms: the 
resonance term 



V2— f — ) cos (2y?i - ^2) 



(73) 



and "swinging" term 



3 1 
T + — 

UJl + 1 UJl 



/1/2 cos {2ipi — 2lp2 



(74) 



the direction m of which is the nearest to the resonance direction r. After that the direct 
application of the destruction stochastic layer criterion leads to the value of critical 
energy in the left well Ecr ~ —0.51 which is in well agreement both with the result of 
numerical simulation and with the predictions of the negative curvature criterion (—5/9) . 
Immediate analysis of integrable part of Hamiltonian (^) shows, that in the right well the 
resonances are absent at negative energies and the transition to large-scale stochasticity 
is possible only with the reach of saddle energy in the complete correspondence with the 
numerical results . 

We can briefly formulate the main results relating to the definition of critical energy 
of the transition from regular motion to stochastic one concluding comparison of the 
results of numerical simulation with the analytical estimates . 

1. The critical energy of the transition to chaos consistents well with the predictions 
of the negative curvature criterion for the PES with unique minimum. 
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2. The critical energy of the transition to chaos either is equal to the minimal energy 
on the zero curvature line or coincides with the saddle energy for the PES with few local 
minima . 

3 The mixed state is observed in the interval f/min {K = 0) < E < Es for the potential 
with several local minima. 

4. The possibility of application of negative curvature criterion (NCC) to the par- 
ticular local minimum is determined by the structure of the Poincare surface of section: 
for any local minimum, 



2.5 BirkhofF-Gustavson normal form 

The structure of the Poincare surfaces of section can be reproduced not resorting to the 
numerical solution of the equations of motion. For this, let us use the method of treating 
non separable classical systems that was originally developed by Birkhoff and later 



was extended by Gustavson [^. The result, obtained by Birkhoff, is the following: if the 
given Hamiltonian H , which can be written as a formal power series without constant 
or linear terms, and such that the quadratic terms can be written as a sum of uncoupled 
harmonic oscillator terms with incommensurable frequencies, then there is a canonical 
transformation that transforms original Hamiltonian into a normal form. The normal 
form is a power series in one- dimensional uncoupled harmonic oscillator Hamiltonian. 
Birkhoff 's method was applied by Gustavson in order to obtain power series expressions 
for isolating integrals, and to predict analytically the Poincare surfaces of section for 
the Henon-Helies system. Since the treated potential had commensurable frequencies, 
Gustavson had to modify Birkhoff's method only somewhat. 

We shall consider the procedure of the transformation to Birkhoff's normal form for 



Hamiltonian which is a power series in coordinates u and momenta v |E7 



H (u, v) = H^^^ (u, v) + H^^^ {u,v) + ... (75) 
where H^^^ has the form of a homogeneous polynomial of degree s: 

//W (m, v)= J2 ciiju'v^s = 2, 3... 

\i\ + \j\=.s 

u' = u{'u^\..u'^ \i\=H+i2 + ... + iN (76) 
For systems in which H^"^^ is positive definite, then there exists a canonical transformation 
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{u,v) {q,p) which transforms H^"^^ into the form 

N 

2 



^ 1 

H^'Hm = J2l-^^{Ql+pl) (77) 

k=l 

Let H{q, p) be Hamiltonian with H^"^^ (g, p) as given in eq. (|77D . Then we say that H{q, p) 
is in a normal form, if 

DH{q,p)=0 (78) 

where 

This is equivalent to the requiring that the Poisson bracket of H^'^^ with H vanish, 
since D is given by 

D = -[H^^\] (80) 

A power series Hamiltonian can be transformed to the normal form by a sequence of 
canonical transformations, where each one reduces the non normalized term of the lowest 
degree to the normal form. The generating function, necessary for the performance of 
individual transformation, is defined by 

k 

The connection between old {p, q) and new (P, Q) canonical variables is 

Qk = qk + ^rj^, Pk = Pk + ^—, H{p,q) = T(P,Q) (82) 
C-nt oqk V / 

where r[P,Q) is the Hamiltonian in the new variables. If we expand H and F in a 
Taylor series about P and q , and then collect and equate all terms of equal degree, the 
following equation for W^^^ is obtained 

DW^'^ (P, q) = r(^) (P, q) - H^'^ (P, q) (83) 
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In order to solve eq. for W^^^ , we make a transformation to variables in which D 
is diagonal: 



Under this transformation D ^ D , where 



_d_ 

dr]k 



It follows that functions of the form 



^l\l2m\m,2 — '/l '/2 SI S2 



are the eigenfunction of D , i.e. 



Iil2mim2 



i^uJk {rrik - k) 



Consequently, 



hl2m-im2 



^ X] {mk - Ik) 



^lil2miT 



Now the equation (|83D may be solved for W^^^ 



(84) 



(85) 



(86) 



(87) 



(89) 



Here H^^^ is a known function. However, so far, T^^^ has been unspecified. Now F*^*^^ 
is determined from the requirement that W^^^ is finite. Clearly, T^^^ must be chosen 
so that to exactly cancel any terms in H^^^ , which would give a vanishing dominator 
in eq. (P5|) . So long as the frequencies are incommensurable, the only terms that must 
appear in T^^^ are those for which = Ik for all k ; those terms are 



(^^i-^i)""' (^^26 



\ 7712 



^ {Pi + ^l) 



(90) 
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Such terms are called null space terms, and the remaining terms are called range space 
terms. Therefore, if H^^^ is separated into null space terms iV(^) and range space terms 

and if we require r*-"^-* to cancel the null space term in H^^^ , eq(|5UD results in 

f (5) ^ ^(S) 

w(s) = n~'R^'^ ^^^^ 

In order to make this solution the unique one, it is sufficient to require the generating 
function containing no null space terms. 

If the transformation to Birkhoff's normal form is made, the harmonic oscillator 
terms can be transformed to action variables via the canonical transformation ([651) . 

We cite as an example the normal form (up to s = 6 ) for the Hamiltonian ( |52D at 
a = 2 in the neighbourhood of right minimum 



H (Ji, h) = 2.613/1 + 2/2 - 0.219/1/2 - 0.017/1^ - 0.375/2^ 
-0.005/1^ - 0.028/1^/2 - 0.122/1/2^ - 0.133/1 



f93) 

3 n noo r n i oo r r2 n i qq r3 v / 



The commensurability of frequencies results into the extension of the set of null space 
terms. As a matter of fact 

and at the realization of condition 

(miUi — 77121^2) = (95) 

D^^rf^^^^'^ would diverge. In order to avoid this occurrence, r^*^-* must be chosen to 
cancel any such additional terms in H'^^'^ . Aside from this change, the procedure is 
identical to that for the incommensurable case. 

The reduction of Hamiltonian to the normal form solves the question about the 
construction of full set of approximate integrals of motion. The latter can be found by 
transformation of the variables of action to initial variables. The solution of equations 

H{Px,Py,x,y) = E 

HPx,Py,x,y) = lo (96) 
X = const 
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allows to find the set of intersections of phase trajectory with selected plane {x = const 
) and by doing so to reconstruct the structure of the Poincare surfaces of section. 

The Poincare surfaces of section for the quadrupole oscillations of nuclei Kr"^"^ , which 
are constructed in such a way are shown in Fig. 8. The Hamiltonian, describing these 
oscillations, up to the terms of the sixth degree with the respect to deformation is the 
following 



The parameters determining the dynamics of the particular nucleus, is calculated for 
the isotopes of Krypton in paper . The Poincare surfaces of section for the equivalent 
quartic potential, the parameters of which are selected from the condition of coincidence 
of the situation of critical points and the values of energy in them are shown in the same 
figure. 

2.6 Regularity-chaos-regularity transition 

The transition regularity-chaos for nonintegrable low-dimensional Hamiltonian, which is 
going on as the energy or the amplitude of external field increases, is a well-investigated 
process. The critical energy of this transition, calculated within the framework of dif- 
ferent scenarios of stochastization, at least for potentials with the simple geometry, is 
in the agreement with the results of numerical simulation. For the systems with local- 
ized region of instability (region of negative Gaussian curvature or region of overlap of 
non-linear resonances) at further increasing of energy, one would expect the return to 
the regular motion; and the critical energy of this new transition chaos-regularity Ecr2 
will be determined by the top boundary of region of instability. For the Hamiltonian of 
quadrupole oscillations (^) 



Recall that for the PES with < < 4 at all energies the motion remains regular. 



^=2 {Pl+Pl)+U{x,y) 




(97) 




(98) 
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Figure 8: The Poincare surfaces of section for different values of energy for the central 
minimum of Kr"^^ . a) obtained by the numerical integration of equations of motion; b) 
obtained with the help of normal forms; c) for the equivalent potential of quadrupole 
oscillations of the fourth degree. 
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Figure 9: Figure represents "the phase diagram" that allows to determine energetic 
intervals of regular and chaotic motions for the fixed value of W . Recall that for the 
PES with < < 4 at all energies the motion remains regular. 

Fig. 9 represents "the phase diagram" that allows to determine energetic intervals of 
regular and chaotic motions for the fixed value of W . 

In the connection with the discussion of possibility of existence of additional transition 
chaos-regularity, let us trace the change of sign and absolute value of Gaussian curvature 
for the PES of Kr'^^ (^) and the equivalent PES of the forth degree (|16]). The parameters 
of the equivalent PES are fitted according to the condition of coincidence of the position of 
extremuma of these two potentials and values of energy in them. The value of Gaussian 
curvature K{x = 0,y) of the PES of Kr'^'^ and the value of the equivalent potential 
of the fourth degree are represented in Fig. 10. We see, that the region of negative 
curvature of the PES of the equivalent potential occupies considerably larger region 
of space and has one order of value larger at y > than for the PES of Kr'^'^ . The 
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Figure 10: The value of Gaussian curvature K [x — Q, y) of the PES of Kr^'^ (sohd hne) 
and the value of the equivalent potential of the fourth degree (dashed hne). The value 
of K in the shaded range of y are very small (less then 0.1). 
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Figure 11: The part of phase space 5*% .with chaotic traectories as a function of the 
energy {Eg - the saddle energy). The solid line - for potential of isotope Kr"^"^ and 
the dashed line - for equivalent potential of the fourth degree. It is seen that for the 
nucleus of Krypton the critical energy of the transition chaos-regularity approximately 
twice exceeds the critical energy for the equivalent potential. With increasing of energy 
in both potentials the 100% chaos appears. However, at energies considerably exceeded 
the saddle energy, and for the nucleus Kr^"^ and for the equivalent potential the regular 
character of motion is restoring. 
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measure of divergence of the classical trajectories, which leads to the rise of the stochastic 
properties in the system, is determined by the size of the region and absolute value of 
negative Gaussian curvature. This circumstance allows us to understand qualitatively 
the reason of transition to chaos in nucleus Kr'''^ at comparatively higher energies, than 
in the equivalent potential: the factors which determine the chaotic character of motion 
in nucleus Kr'^^ are essentially suppressed. Comparatively small region of space, where 
the negative Gaussian curvature of the PES of isotope Kr'^'^ is localized, determines 
the character of motion at energies essentially exceeded the saddle energy. The part of 
chaotic trajectories S, % in relation to all considered trajectories as a function of energy 
for the motion in deformation potential of isotope Kr^'^ and equivalent potential of the 
fourth degree is represented in Fig. 11. It is seen that for the nucleus of Krypton the 
critical energy of the transition chaos-regularity approximately twice exceeds the critical 
energy for the equivalent potential. In both potentials the 100% chaos appears with 
increasing of energy. However, at energies considerably exceeded the saddle energy, and 
for the nucleus Kr'^'^ and for the equivalent potential the regular character of motion is 
restoring. This new transition chaos-regularity is illustrated by the Poincare surfaces of 
section at the super higher energies, which are represented in Fig. 12. 

An earlier appearing of regular motion, at high energies for the quadrupole oscillations 
of nuclear surface of Kr"^^ comparing with the oscillations in equivalent potential, as with 
low energies, is explained by smaller absolute value of Gaussian curvature and larger 
degree of its localization. Underline, that similar reconstruction of regular motion at high 
energies must occur for any potential with localized region of instability. In particular, it 
occurs for isotope of Kr"^^ . For isotopes Kr^^'^^ the negative curvature is not spatially 
localized and that is why the regular character of motion is not reconstructed at energies 
accessible to numerical calculation. 

In conclusion of this section we would like to note the similarity in structure of 
phase space of considered two-dimensional autonomous Hamiltonian system with the 
compact region of negative Gaussian curvature and one-dimensional system with periodic 



perturbation [)48 



The behavior of the width of the resonances, Wk = \ {Wk+i + Wk) , and the distances 
between them, A/^ = \Ik+i ~ Ik\ , as a function of the resonance number is the simplest 
when the satisfaction of resonance overlap condition (2.3.20) for number ki (at a fixed 
level of the external perturbation) guarantees that this condition holds for arbitrary 
k > ki . This is precisely the situation which prevails in the extensively studied systems 
of a ID Coulomb potential |^ and a square well |Q subjected in each case to a 



1/ _2, 

monochromatic perturbation. In the former case we have Wk ~ '6 and A/^. ^ k '3 
while in the latter we have Wk ~ k~^ and Alk ~ [/c (/c + 1)]~^ . 
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Figure 12: The Poincare surfaces of sections illustrating the restoration regular character 
of motion at energies considerably exceeded the saddle energy for the nucleus Kr'^'^ (b,d) 
and for the equivalent potential (a,c). The figures (c,d) at the middle energy: E = lOOEs 
, and (a,b) - at high energy: E — lOOOEs . 
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Figure 13: The resonant spacings AJ^ and the mean widths Wk as functions of the 
resonance numbers k . On the left: for one- dimensional Coulomb, on the right: for 
square-well potential. The critical point ki separates regular range (R) from the chaotic 
one (C) . 



As can be seen from Figure 13 there is a regularity-chaos transition (we will call 
this a "normal" transition) for both ID Coulomb problem and a square well, since there 
exists a unique point ki such that at k > ki the condition Wk > A/^ always holds. The 
motion is therefore chaotic. However, as the behavior of the widths of the resonances 
and of the distances between them as a function of the resonance number becomes more 
complex, we can allow the appearance of an additional intersection point and thus a new 
transition: a chaos- regularity transition, which we will call "anomalous". This is also 
the exotic possibility of the intermittent occurrence of a regular and chaotic regions in 
the phase space. 

We demonstrate that an anomalous chaos-regularity transition occurs in a simple 
Hamiltonian system: an anharmonic oscillator, subjected to a monochromatic perturba- 



tion ||8|. The dynamics of such system is generated by the Hamiltonian 

H {p,x,t) = Ho{p,x) + Fxcosnt (99) 
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with the unperturbed Hamiltonian is 

(p, x) = ^ + Ax" = E {n = 2l,l>l) (100) 

Considered system fills a gap between two extremely important physical models: the 
harmonic oscillator {n = 2 ) and square well (n = oo ). In terms of action-angle 
variables (/,(?) , the Hamiltonian Ho{p,x) becomes 

Ho{i)=(-^iy (101) 



aG (n) 



where 



, , 2V2'KmT + 2n , , 

Ar(i + i) - + 2 

The resonant values of the action 1^ that can be found from the condition ku (Ik) = 
uj{I) = ^ are 

A classical analysis, based on the resonance-overlap criterion, leads to the following 
expression for the critical amplitude of the external perturbation 

per ^ 2(2-3n)/3 1 ^1 ( 9M] f]2n/3^4/3 L(n+2);3 _ ^ ^)(n+2)/3l ' ^^q^^ 



An P Xk \ vr 

where Xk is a Fourier component of the coordinate x{I, q) . Expression (|104|) solves the 
problem of reconstructing the structure of the phase space for arbitrary values of the 
parameters. 

The "phase diagram" in Fig. 14a can be used to determine, at the fixed level of the 
external perturbation, the energy intervals of regular and chaotic motion. The snap-shot 
of E{x) at the right in Fig. 14b confirms that an anomalous chaos-regularity transition 
occurs. We can clearly see isolated nonlinear resonances which persist at large values of 
k , and near which the motion remains regular. The reason for this anomaly is explained 
by Fig. 15. The plots of the resonance widths and of the distances between resonances 
in this Figure demonstrate that there are two intersection points k = ki and k = k2 , 
rather than one. Consequently, there is an anomalous chaos-regularity transition. 
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Figure 14: The "phase diagram" - Figure (a) - the resonance energies Ej^ {k - the reso- 
nance number) against the critical values of the external perturbation for Hamiltonian 
(|10(J|). The upper part of curve (A; = 1 — 5) corresponds to a normal transition regularity- 



chaos but the lower part of curve (A; = 5 — 15) corresponds to an anomalous transition 
chaos-regularity. The snap-shot of E (x) at the right in Figure (b) confirms that an 
anomalous chaos-regularity transition occurs. We can clearly see isolated nonlinear res- 
onances that persist at large values of k , and near which the motion remains regular. 

Thus, one can observe the transition regularity-chaos- regularity just as in the case 
of 2D autonomous Hamiltonian system, so in the case of ID system with periodic per- 
turbation. The reason of the additional transition in both cases is common: localized 
region of instability. In the first case this reason is a localized region of negative Gaus- 
sian curvature, in the second one this reason is a localized region of overlap of non-linear 
resonances. 

2.7 Chaotic regimes in reactions with heavy ions 

Outlined in the previous sections, the common conception, concerning ossible stochas- 
tization of quadrupole nuclear oscillations of high amplitude, is confirmed by the direct 
observations of chaotic regimes at mathematical simulation of reactions with heavy ions 
i- 

The time-dependent Hartree-Fock (TDHF) theory ^7\, ^ constitutes a well-defined 
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Figure 15: The resonant spacings Alk and the mean widths Wk as functions of the 
resonance numbers k for potential (|100|) with n = 8 . One can see two critical points 
ki and /c2 • The second critical point corresponds to an anomalous chaos regularity 
transition. 



starting point for the study of such processes. The TDHF equations can be obtained 
from the variation of the many-body action S , 



t2 



S = dt{^{t) 



.d_ 



H 



(105) 



In this expression H is the many-body Hamiltonian and the A -nucleon wave function 
\E'(t) is chosen to be of determinantal form, constructed from the time-dependent single- 
particle states ipx (f, t) 



^ (fi,f2...fA;t) 



1 

7T\ 



det \\i)x (r, 



(106) 
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The variation of eq.(nJS) is an independent variation with the respect to the single- 
particle states ipx and ■^/'^ and yields the equations of motion, 

and a similar equation for ipl^ (r, t) . 

The classical nature of these equations can be put into a better respective via the 
definition of classical field coordinates 0a if, t) , and conjugate momenta vr^ (r, t) 

/ (108) 

Then the result is the Hamiltonian's equations 



It ~ ^TTA (f,t) (109) 
d-Kxjr^t) _ 5{H) 
dt 5(t)x{r,t) 



The TDHF equation ( |107| ) and its complex conjugate are solved on a three-dimensional 



space time lattice [£2| with initial wave functions of the form 

lim^A (r,t) V2cos (kxr - 6xt^ Xx f 
limTTA {f,t) —>■ v^sin (^\f — e\t^ 



(110) 



where ^a is the solution of the static HF equations hxx (r) = SxXx (r) A = 1, ...A and kx 
is the parameter of the initial boost. 

TDHF calculations for head-on collisions ofHe^ + C^'^, C^^ + C^^O+), and He'^ + Ne^° 
have been performed by Umar et. all. P] at bombarding energies near the Coulomb 
barrier. The results are interpreted in the terms of their classical behavior. The initial 
energy and the separation of the centers of the ions are the parameters labeling the initial 
state. After the initial contact compound nuclear system [O^^ or Mg'^'^ ) relaxes into 
a configuration, undergoing quasi periodic or chaotic motion. The analysis of nuclear 



44 




Figure 16: Poincare phases space plot of M^/ (t) vsM^j (t) for isoscalar quadrupole mode 
for the Mfyf^^ system M . 



density multipole moments <yMLi (t) , Mli (t) j has been apphed for classifying those 
solutions Poincare sections. The definitions of the moments are as follows 



Muit) = / rfVr^FLAf(r)p/(F,t) 

\ (111) 
Mli (uj) = dt exp {-iut) Mu (t) 

where isoscalar (J = ) and isovector (/ = 1 ) densities, 

P'^"^'') { p,{f,t)~p^{f,t) 1 = 1 

The proton pp and neutron p„ densities in terms of the field coordinates (pi and momenta 
TTi , are 



■t) = lYl O^A,, (r, t)\' + \(f)x,, (f, t)\'] ,q=p,n (113) 
The isoscalar quadrupole mode (L = 2,/ = ) is shown in Fig. 16 for the Mg'^'^ 



P,[r,^, 2 

A 
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Figure 17: The autocorrelation function Cli (t) in units Mev i^fm^ sec) as a function 
of time for the isoscalar quadrupole mode in the Mg'^^ system . 

nuclear system seems to be filling most of the available phase space in Poincare section 
jiWi/ (t) ,Mli (^)| • The corresponding autocorrelation function (Fig.l7) 

oo 

/d 
-^^expM)|M2oM|' (114) 
2n 

— oo 

damps fast. All this favours the view that corresponding motion is closer to be stochastic 
rather than quasiperiodic. 
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